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۰ 3 60اه کتمع‌صهل مط) ممتامه هن مه مممتامومز م۳۷ عصنعممهصه و۵ فعمصمتااصا رالی‌تاوهت فمته‌نمی «مصگ 
۲« ما متل صتوحها مط) صا ماه واتعجعل #«مصی ۵۶ »امه عمط ۷۷1 2660 و1 عمصنا عامصط ب«مصو گم عصتاع0مهصه مطا یت۷۲۵۲۵۵۷۵( 
ومحاتله طعنط 2 عصمتاهاه صمصصمسععمصه سمصه . عصتاقعته 240 رومعومص «مصی ‏ صمتاباماتتاونل مط طا ومعصقطم 
۵ 0حه حمتاهتاوه #«عط ج تصواصا ما اطعا بای عتطا رومناووا عطا مصتیملتفمم رم وتعصوجونن فص غانمتگت 
معط صرح (م) مامح رحل-706عع0 مطا ممتصصتعاعل ما و7۵ 2ع1ونرداجر عمج لمح تمامرصته 2 ما عمنامنصطمع) عصنوصهو ع)0ظظ۲۵ 
عط ۳1 عصماح و0مطامصه سعص فصه همتفوهاه عصتمیه فمباله۷ 62012060 وا عصتصتلم رها 0عصمنوه مه صمتاه‌اناصصته مط 
6 81۷۲(۰) 06۱مص محر تامصص«مصوی عمط عم صمتم‌تگلهمن ممتووععع فطل ۵۶ ففباله لقصمجوعو مضه فتاه 
,02121061615 وا صرح [06مصظ فطل ۵۶ عمصهصرموص مط) ماقالده وه تمقتری ما .0عتهوصرمی صه 0مممنالهبه م۵ فمصممتتام 
-2012 204 2011-2012 وتهع۷ ۲م۷2 معط 0۶ د0مننهم ممتاهل‌نله۷ جح صمتاجتطانلجه مطا عم 0عصمگنهجر مت۱۷۵ فممتاهاباصله 
۵ 300۷۷ عط متقصتاوی 0 2560 7۵5 (۷۵01۲10۸1) )مامت بع7مع 0و ۱۷۵1218 عطا مه ,اه تا060و۲۵ ,2013 
6 ۱02۵0۰ 7۵۵ (-رلآنا۳] ۱۱۳۱ و0۳)) عمج متاهنل۲2 عمط ۳0 .ععله۷ ممتاهته: مطا متقصتادع ۵ 240 2762 
۷ 36 ۷7۵5 (2011-2012) 00و ممتاح‌طانای مطا صد ععمصد مصتاه‌آنصته تم 0مطامصط اععها عط عقظ) 7۵0«مطو والتاوعر 
عطا ده تحعمت‌تگلمم مملوعععع1 لهطمعدعو فطل ۵۶ ممیله۷ #۵ فطا فص 0مطامصط صمتاعطاملقه تمامع؟ «2۲66-02ع0 
معط ,(2012-2013) 0مترمم حمتاهلتله۷ عمط ص .4.17 ۵۲ ممصمع0 متام مط فصه 0۰72 ۵۶ ممتافصتصصه‌امل هم احمتم‌تر1ههع 
2 ۷۱ تصمه‌تااهمع ممتووععع۲ اقتاصصه مج مضه 0مطامص تمامع رمل-ععتوع0 امه فطا طاز/۳ 0مطامصه صمتاه‌ااصته مصنار 
۶ عم صاً وااناعع۲ اععها مطا ۲۵۷10۵0 4.38 ۵۶ ممصمتععنل متام مطا فصه 0.51 ۵۶ عمتامصتصماع که احماهم‌تگلههه 

2 260107207 190061 عصا ۵۶ )29909510 
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چکیده 

با توجه به تأثیر مطالعات برف در مدیریت منابع آب و جلوگیری از مخاطرات ناشی از 
بروز سیل و ازآنجاکه مدل‌سازی رواناب ذوب برف به دلیل تغییرات بسیار در 
پراکندگی توده‌های برف با کمبود یا فقدان داده چگالی برف در حوضه آبریز مواجه 
بوده و ایجاد ایستگاه‌های اندازه گیری برف در ارتفاعات کاری سخت و پرهزینه اسنتث: 
در این مطالعه تلاش شده با ابداع معادله‌ای حجدید و با استفاده از تکنیک‌های 
سنجش‌ازدور به روشی ساده‌تر و فیزیکی‌تر پارامتر فاکتور درجه-روز (0) را محاسبه 
کرده و سپس شبیه‌سازی با تعریف مقادیر محاسبه شده آن به روش‌های کلاسیک و 


جدید در کنار مقادیر سالانه و فصلی ضریب فروکش (6) برای مدل رواناب ذوب 


۱ نویسنده مسئول: ۰۹۱۵۱۱۵۶۶۰ 1 
نحوه ارجاع به اين مقاله: 

شهرکی مجاهد. رضا؛ علیزاده امین؛ ثنایی نژاده سید حسین؛ فرید حسینی» علیرضا؛ زرین, آذر؛ ۱۶۰۳. مدل‌سازی رواناب ذوب 
برف توسط مدل ]۹1۹1۷ و برآورد پارامتر فاکتور درجه-روز با استفاده از تصاویر ماهواره‌ای تابش خالص (مطالعه موردی: حوضه 


آبریز کارده). جغرافیا و مخاطرات محیطی» ۱(۱۳). صص ۱-۲۲ 
0 و ( 


۲ جغرافیا و مخاطرات محیطی شمارة اول 


برف (51801) انجام شده و نتایج حاصل مورد ارزیابی و مقایسه قرار گرفته است. به- 

منظور ارزیابی عملکرد مدل و پارامترهای آن. شبیه‌سازی برای دوره‌های واسنجی و 

درستی‌سنجی به‌ترتیب برای سال‌های آبی ۲۰۱۱۲۰۱۲ و ۲۰۱۲۲۰۱۳ انجام شد و از 

محصول پوشش برف (0010۸1) ۷1015 جهت برآورد سطح پوشش برف و 

برای برآورد مقادیر تابش خالص از محصول تابش خالص -056۴5) ۱۴۵0 

(۱۳1۳[/1-۲ استفاده گردید. نتایج نشان داد بهترین روش شبیه‌سازی رواناب در 

دوره واسنجی (۲۰۱۲-۲۰۱۱) استفاده از روش محاسبه فاکتور درجه-روز جدید و 

کاربرد دو مقدار ضریب فروکش جریان فصلی با ضریب تبیین ۰/۷۲ و درصد اختلاف 

خی ۶/۱۷ ای جر فوره اففازی 0۲۳۱۳۵۲۸۲ ند شبه‌سارش: وزاتاب: نا 

روش محاسبه فاکتور درجه-روز جدید و یک ضریب فروکش سالانه با ضریب تبیین 

۱ و درصد اختلاف حجمی 1/۳۸ بهترین نتیجه را از نظر ارزیابی معیارهای دقت 

مدل ارائه کرد. 

کلیدواژه‌ها: ذوب برف. فاکتور درجه-روز» سنجش‌ازدور 1۷10115 تابس خالص 
۱- مقدمه 

ایران سرزمین کم‌آبی است که در منطقه خشک و نیمه‌عشک کره زمین قرار گرفته است و لذا بارش برف 
مهم‌ترین عامل تأثیرگذار در میزان ذخایر آب به‌ویژه در فصول گرم سال محسوب می‌گردد (قنبرپور و 
همکاران ۱۳۸۶). برف یک شکل از بارش است که به دلیل تأخیر در تبدیل آن به رواناب تفاوت ویژه‌ای با 
دیگر اجزاء دارد (عالی‌نژاد و همکاران» ۲۰۱۳). رواناب حاصل از ذوب برف به دلیل نقش تأخیری آن منبع 
اصلی تغذیه سفره‌های آب زیرزمینی و در برخی از موارد به دلیل هم‌زمانی با بارش‌های بهاره منشأً بروز 
سیلاب‌های مخرب با حجم جریان بیش از ظرفیت رودخانه‌ها می‌شود (میریعقوب‌زاده و قنبرپور ۱۳۸۹). 
برآورد دقیق آب معادل برف و همچنین تغییرات سطح پوشش آن. در بخش‌های کشاورزی, انرژی؛ 

مدیریت مخزن و هشدار وقوع سیل مژثر است (جهان‌بخش اصل و همکاران. ۱۳۹۵). شبیه‌سازی رواناب 
ذوب برف در مناطق کوهستانی نه‌تنها برای ارزیابی مخاطرات سیل ناشی از ذوب برف در بهار و تابستان 
بلکه برای توسعه و مدیریت منابع آب در حوضه اهمیت زیادی دارد (هاو و همکاران» ٩۲۰۱؛‏ اما در 
کشورهای درحال‌توسعه و مناطق دورافتاده با کوه‌های مرتفع» برآورد برف روی سطح زمین به دلیل 


3 2 ۰ 1 
توپوگرافی پیچیده عدم دسترسی و محدودیت‌های مالی بسیار سخت و محدود است (ویوویچ و جاکوبز ‏ 


1 ۲۱۸0 6 1. 
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سال سیزدهم مدل‌سازی رواناب ذوب برف توسط مدل ٩16۷۲‏ و برآورد پارامتر فاکتور درجه-روز با استفاده از... ۳ 


۱+ با وجود تغییرات زیاد در توزیع زمانی و مکانی برف انباشته شده. اطلاعات مشاهده‌ای زمینی برف 
اندک و معمولا به نواحی با ارتفاع کمتر محدود می‌شود که این موضوع کاربرد سنجش‌ازدور در تأمین بخشی 
از اطلاعات برف را اجتناب‌ناپذیر نموده است. بعلاوه مدل‌های پیش‌بینی رواناب حاصل از ذوب برف نیاز به 
اطلاعاتی دارند که اين اطلاعات با فقر آمار در مناطق برف‌گیر حوضه‌های آبریز عمدتاً قابل‌دسترس و یا 
تأمین نیستند که با این توصیف تکنیک سنجش ازدور به‌عنوان ابزاری جدید و مفید برای به‌دست آوردن داده- 
های برف مطرح می‌باشد (نجف‌زاده و همکاران. ۱۳۸۳). اندازه‌گیری‌های میدانی که به‌ویژه در مناطق 
کوهستانی انجام می‌شود اغلب بدون استفاده از فناوری‌های جدید غیرممکن است (شهابی و همکاران. 
۶ فناوری سنجش‌ازدور با آماده‌سازی لایه‌های اطلاعاتی پیوسته و با دقتی بالاتر و هزینه کمتر می‌تواند 
خلاً آمارهای هیدرولوژیکی را پر کند (طالبی اسفندرانی و همکاران. ۱۳۹۰). به‌طور کلی. مدل‌سازی ذوب 
برف و رواناب حاصل از آن با استفاده از دو روش انجام می‌شود: روش تعادل انرژی و روش درجه-روز دما 
که روش دوم می‌تواند در غیاب داده‌های کافی مفیدتر باشد (سنزبا و همکاران؛ ۲۰۱۵). پارامترهای برف نقش 
بسیار مهمی در چرخه هیدرولوژیکی, الگوی آب‌وهوای محلی و تغییر اقلیم ایفا می‌کنند و درنهایت بر کل 
رواناب جریان رودخانه‌ها و آب در دسترس از طریق مقدار و زمان ذوب برف تأثیر می‌گذارند (نانچهانی و 
همکازان ۲۲۳ 

در مرور پژوهش‌های مرتبط انجام شده در این حیطه. مارتینک (۱۹۸۹) در سایت آزمایشگاهی 
طاممزدب1 ۷1551 واقع در الا سوئیس میزان ذوب ساعتی برف را در سراسر یک فصل ذوب برف محاسبه 
کرد و بده پوشش برف به‌طور مداوم توسط یک لایسیمتر برف اندازه‌گیری شد. مقادیر به طی تمام دوره 
فرسایش بیشتر در محدوده ۰/۲ تا ۰/۲۵ به‌دست آمد و مقادیر محاسبه شده می‌تواند برای مطالعات حرکت 
آب در پوشش برف استفاده شود. کوستاس و رانگو (۱۹۹۶) یک مولفه تابش ساده را با روش درجه-روز در 
تلاشی برای بهبرد ی آوزد خوی بیرف یل ۵9 کیب گروون هلگ از انم مطالته: انساد قاکون ذوت 
برف فیزیکی‌تر برای 515۷1 توسط ترکیب روش بودجه تابش ساده شده با روش درجه-روز بوده است. برای 
روش درجه-روز محدود مقدار ثابت ۰/۲ برای 0 در طول فصل ذوب برف استفاده شد؛ در حالی که در 
روش درجه-روز اصلی مقادیر 0 به‌تدریج از ۰/۶۸ در ماه مه افزایش یافت و به ۰/۵ در ژوئن و ۰/۵۲ در ماه 
ژوئیه رسید. بروبکر و همکاران" (۱۹۹7) پژوهشی در حوضه آبریز اسلیپرز ریور انجام دادند. هدف از اين 


مطالعه ت کیب ورودی تابش در مدل روائاب ذوب برف بود. مقدار فاکتور درجه-روز محدود شده (0) برای 
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تمامی ماه‌های فصل ذوب در تمام سال‌های مورد مطالعه مقدار ۰/۲ برآورد شد. نتایج به‌دست آمده نشان داد 
که نسخه ترکیبی تابش عملکرد بهتری نسبت به نسخه درجه-روز متداول 8181 داشته است. وفاخواه و 
همکاران (۲۰۱۶) مدل تابشی و درجه-روز 58۷ را مقایسه کردند. آن‌ها در مدل تابشی 101 از فاکتور 
درجه-روز محدود شده (,0) استفاده کردند و در این مطالعه. تابش خالص با استفاده از تصاویر ۷10۳16 
برآورد شد. در نتایج به‌دست آمده مشخص شد ضریب فروکش و ضریب رواناب ذوب برف حساس‌ترین 
پارامترها هستند و مدل تابشی 88۷1 نتیجه بهتری با توجه به معیارهای دقت به‌دست آورد. استیلی و 
همکاران" (۲۰۱۷) محصولات برفی ۷0116 برای مدل‌سازی رواناب ذوب برف را در حوضه آبریز رودخانه 
ریوگرانده مورد مطالعه قرار دادند. آن‌ها به کاربران توصیه می‌کنند هنگام شبیه‌سازی ذوب برف با مدل 518۷ 
یکی از محصولات اصلاح شده با سایبان "180۳50۸0 یا محصولات مخصوص به موقعیتشان را از 
مجموعه 0010۸1 ]5 6 001166102 انتخاب و استفاده کنند. شی و همکاران" (۲۰۱۸) در مطالعه‌ای 
الگوریتم‌های بهینه‌سازی را برای کالیبراسیون پارامترهای متغیر در زمان مدل رواناب برف (518۷) در شبیه- 
سازی جریان آب ناشی از ذوب برف برای حوضه رودخانه ماناسی بررسی کردند و نشان دادند برای شبیه- 
سازی یک فصل ذوب برف. طول زیردوره تأثیر قابل‌توجهی بر عملکرد مدل دارد. نانج‌های و همکاران 
(۲۰۲۰) در پژوهشی تغییرات مکانی و زمانی پارامترهای برف در حوضه دیبانگ» آروناچال پرادش, هند را با 
استفاده از مدل رواناب ذوب برف توزیع شده فضایی (5511) همراه با سنجش‌ازدور و تکنیک‌های 615 
را مورد تجزیه‌وتحلیل قرار دادند و نتیجه گرفتند که می‌توان آن را برای استفاده در حوضه‌های اندازه‌گیری 
نشده یا کمی اندازه‌گیری شده در منطقه هیمالیا توصیه کرد. زیمنس و همکارانش" (۲۰۲۱) تخلیه رواناب در 
حوضه بالادست رودخانه آتاباسکای در غرب کانادا را با استفاده از مدل رواناب ذوب برف ]812 و 
مشاهدات سنجش ازدور برف ۷10115 را شبیه‌سازی کردند. نتایج نشان داد گنجاندن تغییرات مکانی و زمانی 
در فاکتور درجه-روز و جداسازی آن برای مناطق یخ‌زده و بدون توده برف برای بهبود شبیه‌سازی شرایط 
مختلف برف در طول چند سال مهم بوده است. همچنین آنها اظهار داشتند این مدل در شبیه‌سازی جریانات 
تاریخی به‌خوبی عمل کرده است. 

میریعقوب زاده و همکاران (۱۳۹۰) در حوضه آبخیز سد کرج نشان دادند که مدل 911 با ضریب تبیین 
۷ و ۰/۹۶ به‌خوبی قادر به مدل کردن فرآیند رواناب حاصل از ذوب برف می‌باشد. در پژوهشی دیگر 
ابراهیمی و همکاران (۱۳۹۵) با استفاده از مدل 5121 و تصاویر هشت روزه سنجنده 1۷۲01715 رواناب روزانه 
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سال سیزدهم مدل‌سازی رواناب ذوب برف توسط مدل ٩16۷‏ و برآورد پارامتر فاکتور درجه-روز با استفاده از... ۵ 


حاصل از ذوب برف را برای حوضه آبریز نهاوند. مورد شبیه‌سازی و ارزیابی قرار دادند. بالاترین ضرایب 
همبستگی و تفاضل حجمی برای مدل در دوره واسنجی به‌ترتیب ۰/۷۵ و ۳/۱۲- و در دوره درستی‌سنجی 
پرابر ۰/۷۹ و ۲۰/۳۷ به‌دست آمد. تحلیل حساسیت پارامترهای ضریب درجه-روز. ضریب رواناب برف و 
باران مدل انجام گردید که نشان دهنده حساسیت کمتر مدل به افزایش و کاهش ضریب رواناب باران در میان 
سایر پارامترها بود. اربابی سبزواری و نجاری (۱۳۹7) به پژوهشی مشابه در ارتفاعات حوضه آبخیز طالقان 
پرداختند. مدل» آبدهی روزانه را با مجذور ضریب همبستگی 40 شبیه‌سازی کرد. نتایج این تحقیقات نشان 
می‌دهند که مدل شبیه‌ساز 5181 ابزاری بسیار مناسب برای برآورد میزان رواناب حاصل از ذوب برف می- 
باشد و از دقت بالایی برخوردار است و این مدل جهت مدیریت منابع آب؛ کارامد. غي‌باشد. دز خوهشی 
مشابه. رشیدی و همکاران (۱۳۹۲) با استفاده از مدل 83۷ رواناب روزانه حاصل از ذوب برف را با 
استفاده از تصاویر سنجنده ۷0115 در حوضه آبریز دربند سملقان محاسبه کردند. نتایج کاربرد مدل شبیه- 
سازی موفق و قابل قبولی را نشان داده است. نتایج برای تحلیل حساسیت مدل به پارامترهای ضریب رواناب 
برف. فاکتور درجه-روز و افت محیطی دما نشان داد که سه پارامتر یاد شده در بالا بر ضریب تبیین و درصد 
تفاضل حجمی مژثرند که البته مدل نسبت به تغییرات افت محیطی دما نسبت به دو پارامتر دیگر از حساسیت 
کمتری برخوردار است. در پژوهشی به‌منظور پایش زمانی و مکانی سطح پوشش برف. میزان ضریب تبیین 
(87) به‌دست آمده برای مدل‌های ]510 و "713۷ در مرحله واسنجی به‌ترتیب برابر ۰/۷۱ و ۰/۱۱ و مقدار 
۲ و ۰/۹ برای مرحله اعتبارسنجی به‌دست آمد (ارتیمانی و همکاران ۱۳۹۸). مدل 512۷۲ به دلیل استفاده 
از تصاویر ماهواره‌ای» نسبت به مدل 1113۷ عملکرد قابل قبول‌تری در شبیه‌سازی رواناب حاصل از ذوب 
برف از خود نشان داد. عالی‌نژاد و همکاران (۱۶۰۰) رواناب ناشی از ذوب برف را در حوضه‌ی قره‌سو با 
مدل ]۹1 شبیه‌سازی کردند. نتایج شبیه‌سازی بیانگر دقت بالای این مدل می‌باشد. به‌طوری‌که ضریب تبیین 
(2) برای سال‌های آبی ۹۵-۹7 الی ٩-۹۷‏ به‌ترتیب معادل ۰/۹۳ و ۰/٩‏ و درصد حطای حجمی (,2) آن نیز 
به‌ترتیب ۰/۳ و ۳/۳۳ به‌دست آمد. 

فقر آمار و اطلاعات در حوضه‌های آبریز یکی از چالش‌های جدی پیش روی مدل‌سازی رواناب ناشی از 
ذوب برف می‌باشد که کار تعیین و محاسبه پارامترهای ورودی به مدل را سخت و دشوار می‌سازد. یکی از 
پارامترهای ورودی مهم در مدل رواناب ذوب برف (5[8) فاکتور درجه-روز (0) است. مقادیر فاکتور 
درجه-روز (0) در حوضه‌های آبریز به‌صورت متداول و کلاسیک توسط داده‌های اندازه‌گیری شده چگالی 
برف محاسبه و تعیین می‌شود و سپس به مدل 51۷۲ وارد می‌گردد. اين داده‌ها به دلیل توزیع مکانی و زمانی 
سطوح پوشیده از برف و کمبود ایستگاه‌های اندازه‌گیری در ارتفاعات مختلف دارای کیفیت پایین» کمبود و 
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فقدان در زمان و مکان‌های مختلف در حوضه‌های آبریز است. ازآنجاکه برف در مناطق صعب‌العبور تجمع 
می‌یابد. ایجاد ایستگاه‌های اندازه‌گیری و برف‌سنجی در ارتفاعات کاری سخت. زمان‌بر و پرهزینه است و با 
توجه به اینکه در مدل 518۷ عمق ذوب برف روزانه رابطه مستقیم با پارامتر فاکتور درجه-روز و دما دارد و 
از سوی دیگر ملفه تابش خالص در محاسبات فاکتور درجه-روز (0) به روش متداول و کلاسیک آن در نظر 
گرفته نمی‌شود و در نتیجه تأثیر آن در عمق ذوب برف منظور نمی‌گردد. هدف از این مطالعه آن است تا با 
استفاده از معادله‌ای جدید و برگرفته از معادلات گذشته تعیین عمق ذوب برف ارائه شده توسط کوستاس 
(۱۹۹۶) و بروبکر (۱۹۹1) همراه با استفاده از تصاویر ماهواره‌ای تابش خالص و تکنیک‌ها سنجش‌ازدور به 
روشی ساده. کم‌هزینه دقیق‌تر و فیزیکی‌تر فاکتور درجه-روز (0) را محاسبه کرده و سپس مقادیر به‌دست 
آمده را جهت شبیه‌سازی رواناب ذوب برف برای مدل 88۷1 تعریف کرده و عملکرد روش جدید را در 
مقابل روش متداول و کلاسیک محاسبه فاکتور درجه-روز (0). در کنار لحاظ کردن مقادیر متفاوت پارامتر 
ضریب فروکش (6) در شبیه‌سازی رواناب ذوب برف توسط مدل را ارزیابی و بررسی کنیم. در این مطالعه 
نتایج حاصل از شبیه‌سازی‌های مختلف رواناب ذوب برف توسط مدل 518۷ با تعریف مقادیر فاکتور درجه- 
روز (0) محاسبه شده با روش جدید و کلاسیک و تعریف مقادیر سالانه و فصلی برای پارامتر ضریب 
فروکش (6) مورد ارزیابی و مقایسه قرار گرفته است و بهترین روش محاسبه فاکتور درجه-روز (0) و 
بهترین شیوه شبیه‌سازی رواناب ذوب برف در حوضه آبریز مورد مطالعه توسط مدل 51 تعیین گردیده 


است. 
۲- منطقه مورد مطالعه 


حوضه آبریز سد کارده با مساحتی حدود ۵۶6/۲۱۳ کیلومتر مربع در شمال شرق ایران و شمال استان 
خراسان رضوی و در شمال شهرستان مشهد قرار گرفته است. موقعیت جغرافیایی آن بین "۲۵ ۵۹٩‏ تا "۶۵ 
۳ طول جغرافیایی و ۳۸ ۳۹ تا ۵۵۳ ۵۷ ۳۷ عرض جغرافیایی قرار گرفته و این حوضه از زیر حوضه- 
های اصلی حوضه آبریز کشف‌رود به‌شمار می‌رود. کم ارتفاع‌ترین نقطه حوضه در محل سد کارده با ارتفاع 
۰ متر و بلندترین نقطه آن در ارتفاعات شمال خرکت با ارتفاع ۲۹۷۷ متر واقع شده است. 
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شکل ۱- موقعیت جغرافیایی منطقه و ایستگاه‌های مورد استفاده در حوضه آبریز کارده 


۳- مواد و روش‌ها 


مدل روانات ذوب برف (۷7[) که توسط مارتینک در تیال ۱۹۷۵ توسعه داده سل نگ روش مفهومی 


نیمه‌توزیعی درجه-روز دماه برای شبیه‌سازی و پیش‌بینی رواناب روزانه حاصل از ذوب برف و بارندگی در 
تفای کرسیدایی ان مشک و رای ما اناهای 0۱ کاس ربب تشگنا کی 


درجه-روز و ضرایب رواناب برای برف و باران می‌باشند که می‌توانند شبیه‌سازی رواناب را اصلاح و تعدیل 


1 (۷2۲010606 ۰ 
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کرده و بر عملکرد مدل تأثیر بگذارند (پاندی و همکاران؛ ۲۰۱۶). 518۷1 در گام زمانی روزانه و شبیه‌سازی 
جریان مبتنی بر تغیبرات دمای روزانه. بارش و پارامترهای مختلف به نمایندگی از شرایط فیزیکی حوضه 
عمل می‌کنند (الیاس و همکاران؛ ۵ فراهم کردن داده‌های عمده برای شبیه‌سازی و مدل‌سازی جریان در 
حوضه‌های کوهستانی و با توجه به نقش مهم رواناب ذوب برف و همچنین کمبود داده در اين مکان‌ها یک 
چالش است (فیروزی و صادقیان. ۲۰۱۳). همچنین هرگونه تقویت در اساس فیزیکی روش درجه-روز استفاده 
از آن را در پیش‌بینی رواناب ذوب برف و ارزیابی اثرات تغییر اقلیم بهبود خواهد بخشید (کوستاس و رانگو" 
۹ 

مدل رواناب ذوب برف (5[60۷1): رواناب حاصل از ذوب برف و باران را به‌طور روزانه محاسبه کرده و با 
لحاظ کردن فروکش جریان. دبی خروجی از حوضه را برآورد می‌کند. مساحت و محدوده حوضه و ناحیه- 
های ارتفاعی بوسیله نقشه "9۳ منطقه و توسط نرم افزار ۸86615 به‌دست آمد. منطقه به سه ناحیه 


ارتفاعی با فواصل ۰۰۰ متر تقسیم شد (جدول ۱). 


حدول ۱- مشخصات ناحیه‌های ارتفاعی تعیین شده در حوضه آبریز کارده 


مساحت طبقات ار تفاعی ناحیه ارتفاعی 
,۹ ۹9 ( 
۳/۹ ۱-۹۲( ۱۷۰۰۰ 
٩0۸/۹۵‏ ۳۳/۸۲ ۲۳۰۰-۸۰ 
۱1/۳ ۷1۸/۲ ۲۹۰۰-۲۰ 
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شکل ۲- نقشه طبقه‌بندی ارتفاعی حوضه آبریز کارده 


متوسط بارش روزانه از ایستگاه‌های هیدرومتری و باران‌سنجی موجود در منطقه برای سال‌های 
۲ ۱۱ و ۲۰۱۳ ۲۰۱۲ استفاده شد (جدول ۲) و توسط گرادیان بارش برای متوسط ارتفاع 
هییسومتریک هر ناحیه به‌دست آمد. دما متغیر مهمی در محاسبه عمق ذوب برف روزانه است. با توجه به 
پيشنهاد توسعه دهندگان مدل مبنی بر اينکه یک ایستگاه دمایی خوب (حتی اگر در خارج از حوضه قرار 
داشته باشد) به چندین ایستگاه کمتر قابل اطمینان ترجیح داده می‌شود. به این منظور از داده‌های دمای بيشینه 
کته ایستکاه‌های سینوپتیک گلمکان و مشهد استفاده شد (جدول ۳) و توسط گرادیان دمایی برای متوسط 


ارتفاع هیپسومتریک هر یک از ناحیه‌های ارتفاعی استخراج شد. 


۱۰ جغرافیا و مخاطرات محیطی شمارةٌ اول 


حدول ۲- مشخصات ایستگاه‌های هیدرومتری و باران‌سنجی مورد استفاده در حوضه آبریز کارده برای سال‌های 
۱۲ و ۲۰۱۳ ۲۰۱۲ 


ارتفاع از سطح دریا طول حغرافیایی عرض جغرافیایی نوع ایستگاه نام ایستگاه 
9 
٩۵۳ ۱۳۳۲‏ ۳۳۹ هیدرومتری کارده بالادست 
۱۳۷۹ ۹9۹8 ۳0۳۷۲ باران‌سنجی سد کارده 
۱1۹ ۲ 9۵ ۳ 9۳ باران‌سنجی آل 
٩۵4 ۲ ۱۳۳‏ هیدرومتری و باران‌سنجی جنگ 
٩۵۹۳۳ ۱91۹‏ و باران‌سنجی گوش‌بالا 
۱۸۳۰ 0 0 باران‌سنجی مارشک 
٩۵4 ۲ ۱۹۶۱‏ و باران‌سنجی بلغور 


حدول ۳- مشخصات ایستگاه‌های هواشناسی مورد استفاده در حو ضه آیریة کارده برای سال‌های ۳ و 


۳ ۱۲۰۱ 
ارتفاع از سطح دریا (۳0) طول جغرافیایی عرض جغرافیایی نوع ایستگاه نام ایستگاه 
۹۹۹ ۳۸ ۹۵ رت سینوپتیک مشهد 
٩۵۹۷۳ ۱۳۹‏ ۳۹۳ سینوپتیک گلمکان 


سطح پوشش برف یکی از متغیرهای اصلی ورودی به مدل 518۷1 است که به‌صورت روزانه مورد نیاز 
مدل می‌باشد. برای استخراج سطح پوشیده از برف برای سال ۲۰۱۱۲۰۱۲ جهت واسنجی مدل و سال 
۳ برای درستی‌سنجی مدل از محصول پوشش برف روزانه سنجنده ۷01۲5 از ماهواره 16772 
به‌نام ۷10110۸1 استفاده گردید. در این محصولات برای استخراج سطح پوشش برف از الگوریتم 90017 
مهد استفاده شده که در این الگوریتم از شاخص تفاضل نرمال شده پوشش برف ( ۷۳51 استفاده شده 
است. تصاویر به‌منظور تبدیل سیستم مختصات. عملیات برش تصاویر برای هریک از سه ناحیه ارتفاعی 
تعیین شده در حوضه آبریز کارده و محاسبه اطلاعات آماری به محیط نرم افزار ۳۳۷۲ فراخوانده شد و 


۵0۱۷ ۲1]]6۲۵96۵ 0عوتلمصصورم 1 


سال سیزدهم مدل‌سازی رواناب ذوب برف توسط مدل ٩16۷‏ و برآورد پارامتر فاکتور درجه-روز با استفاده از... ۱۱ 


1 
8 مر 
96ص 
4 (وء 
وه 3" 
0 . 
شم هی ی و یش 
لا لام لا سم 
کی کی کي کم ۳ 
۱ ار ۳ 
ناحیه ۸ کچ و ما و مور ی 
ناحیه 8 و زمان (روز) 
ناحیه » تک وک بت 


شکل ۳- (الف) تغییرات سطح پوشش برف در سال ۱۲ و (ب) تغییرات سطح پوشش برف در سال 


۲۰۱۲.۳ در حوضه آبریز کارده 


ضریب رواناب برای کل حوضه آبریز محاسبه و در طی مرحله واسنجی مقادیر بهینه آن به‌دست آمد. 
دمای بحرانی (ج1) نیز طبق توصیه توسعه‌دهندگان مدل, دمای ۲/۵ درجه سلسیوس برای ماه‌های سرد 
سال و دمای ۰/۷۵ درجه سلسیوس برای ماه‌های گرم سال در نظر گرفته شد. گرادیان دما (۲) مقدار پیش 
فرض مدل و مقدار ۰/0۵ درجه سلسیوس به ازای هر ۱۰۰ متر ارتفاع و با توجه به روش انتخابی زمان تمرکز 
در مطالعات فیزیوگرافی که روش کرپیچ می‌باشده زمان تأخیر (1) برابر 7/۵" محاسبه گردید. 

فاکتور درجه-روز (0): فاکتور درجه-روز [17 0۳" طع]» تعداد درجه-روز [1].0 را به عمق ذوب 
برف روزانه [۷1]0 تبدیل می‌کند. فاکتور درجه-روز (0) در خارج از مدل محاسبه شده و سپس مقادیر آن 
برای مدل 8181 تعریف می‌گردد. فاکتور درجه-روز با در دسترس بودن مقادیر اندازه‌گیری شده چگالی 
برف» به شیوه کلاسیک و متداول از رابطه تجربی زیر به‌دست می‌آید. 

9-1. (۲ 

که در اين رابطه ,۵ چگالی برف و 0 چگالی آب می‌باشد. فاکتور ذوب به‌صورت مقدار عمق ذوب شده 

به تعداد درجه-روز بوده و به‌شکل رابطه (۲) بیان می‌گردد. در تلاشی بروبکر و همکاران (1۹۹7) تابش را در 

مدل ترکیب کردند. در این روش معادله ذوب برف به‌صورت رابطه (۳) بیان می‌شود. در این پژوهش با قرار 

دادن معادلات شماره (۲) و (۳) در یک دستگاه معادله جدیدی جهت برآورد فاکتور درجه-روز (0) به‌دست 
آمده است. معادله فاکتور درجه-روز (0) جدید به‌صورت رابطه شماره (۶) می‌باشد: 


۲( 0.1 تب ۳ 
۳( ۳ لا ِ] 11۳۳0۰ تب 1 


زگ 
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۳ ۱10۰ ۳0 


1 0 
1۵ ۱ 


0 


که در این معادله 100 ثابت فیزیکی تبدیل انرژی به حجم یا عمق آب است (4۵۷ ۳/۷۷ 0 ۱/۰۲۱ به 
ضریب درجه-روز محدود شده است (02۷ /610). و" شاخحص درجه حرارت است. 5 تابش خالص 
روزانه است. به دلیل اینکه تغییرات فاکتور درجه-روز (0) به‌صورت ۱۵ تا ۳۰ روزه توصیه می‌شود برای 
محاسبه تابش خالص از محصول تابشی خالص ۸ روزه (۱0۳21۳1/6-5- 0۳8۳5)) 11520 با قدرت تفکیک 
مکانی ۰/۲۵ درجه استفاده شد. جهت محاسبه تابش خالص تصاویر در نرم افزار ۳(۷1 فراخوانده شد و 
مقادیر آن استخراج گردید. با توجه به نتایج مطالعات مارتینک در سال ۱1۹۸۹ کوستاس و رانگو (1۹۹4) و 
بروبکر و همکاران (۱۹۹7) جهت برآورد مقدار ضریب درجه-روز محدود شده (,0) در اين مطالعه نیز مقدار 


۰/۲۵ برای این ضریب در محاسبات مورد استفاده قرار گرفت. 


جدول ۶- مقادیر پارامتر فاکتور درجه-روز (0) محاسبه شده با معادلات کلاسیک و جدید در حوضه کارده 


مقدار واقعی فرحه محدوده نرمال نماد نام متغیر 

7۱-۹ محاسبه فاکتور (0) با معادله کلاسیک 
1 ۱ 2 

۶ - ۰/۲۵ ۳ ۱ (می»/ 3 محاسبه فاکتور (0) با معادله جدید 


ضریب فروکش یک پارامتر بسیار مهم در مدل 8181۷ است. تحلیل داده‌های سری زمانی دبی اندازه‌گیری 
شده روش خوبی برای تعیین 16 می‌باشد باید توجه داشت که 16 ثابت نیست و با کاهش دبی مقدار آن 
افزايش می‌یابد. مدل 516۷1 خود 1 را به‌عنوان ورودی نمی‌پذیرد و بایستی ضرایب ‏ و ل به مدل داده شود 
به این منظور مقادیر دبی جریان در هر روز ببم در مقابل دبی روز قبل :0۵ در قالب محور لگاریتمی رسم 


شد و با استفاده از تموذار هکت آمد+ مقادیر ضرایب ‏ و 7 محاسبه شد. 


حدول ۵- مقادیر ضرایب * و ۷ به‌دست آمده برای معادله فرو کش حریان (16) 


۷ 4 
۱/۷ ۰*۸۹ ماه‌های سرد 
۹۷« ۱/۸۹۵۲ ماه‌های گرم 
۰/۰۹ ۱ سالیانه 


ک20 0056۲۷ طاتقظ ممول [ 
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تا سیزدهم مدل‌سازی رواناب ذوب برف توسط مدل 7 و برآورد پارامتر فاکتور درجه-روز با استفاده اتمه ۱۳ 


مدل ]۹180۷ هیدروگراف رواناب اندازه‌گیری شده و محاسبه شده را به‌صورت گرافیکی نمایش می‌دهد. در 
عین حال برای ارزیابی عملی‌تر از چگونگی انجام شبیه‌سازی از دو معیار دقت. ضریب تبیین (8) و درصد 
تفاضل حجمی " (,) استفاده می‌کند. مقادیر تعیین شده پارامترهای مورد نیاز مدل تا حد امکان براساس 
داده‌های موجود برآورد گردید و در صورت نبود اطلاعات و آمار کافی از مقادیر پیشنهادی توسعه‌دهندگان 
مدل و مطالعات مشابه استفاده و در مرحله واسنجی مقادیر بهینه آن‌ها تعیین شد. در این پژوهش مدل رواناب 
ذوب برف (*51) برای مقادیر به‌دست آمده فاکتور درجه-روز (0) از دو روش کلاسیک (معادله شماره ۱) 
و روش جدید (معادله شماره ۶) در کنار مقادیر متفاوت ضریب فروکش جریان () به‌صورت یک مقدار 
ثابت برای تمام طول دوره شبیه‌سازی و استفاده از دو مقدار فصلی برای ماه‌های گرم و سرد اجرا شد و نتایج 
آن‌ها برای تعیین بهترین روش شبیه‌سازی و عملکرد بهینه مدل توسط معیارهای ارزیابی دقت مدل مورد 


6- نتایج و بحث 


واسنجی مدل: مطابق آنچه شرح داده شد اطلاعات و داده‌های مورد نیاز برای اجرای مدل رواناب ذوب 
برف ۷۷1051 تعریف گردید و سپس مدل برای ماه‌هایی از سال که پوشش برف در سطح حوضه وجود 
ره براغسال اب ۳۱۲ ۳ سصهی واتی اضرا مت ول انسعات سای ۳۱۹۱۱۰۲۰۱۲ شام شش 
مدل» متوسط بارش برف ۱/۷ سانتی‌متر در این سال است که نسبت به دهه قبل و بعد از آن بیشتر بوده و 
باعث عملکرد بهتر مدل در شبیه‌سازی رواناب ذوب برف می‌شود. شکل ۶. هیدروگراف‌های شبیه‌سازی شده 
و اندازه‌گیری شده روائاب ذوب برف را در سال ۲۰۱۱-۲۰۱۲ که توسط مدل 51808 ترسیم شده را نشان 
می‌دهد. که برای ارزیابی و مقایسه بصری بهتر هیدروگراف‌های روش‌های مختلف در کنار یکدیگر آورده 
شده‌اند. براساس معیارهای ارزیابی دقت مدل (جدول 1). اجرای مدل رواناب ذوب برف با فاکتور درجه- 
روز جدید و استفاده از دو ضریب فروکش فصلی بهترین عملکرد را در شبیه‌سازی رواناب در طول دوره 
شبیه‌سازی داشت. در این اجرا. ضریب تبیین (*8) برابر با ۰/۷۲ و درصد اختلاف حجمی (,12) برابر با ۶/۱۷ 
بود. میزان حجم رواناب اندازه‌گیری شده ۷/۵۵۰ میلیون مترمکعب و متوسط رواناب اندازه‌گیری شده ۰/6۸ 
مترمکعب در ثانیه بود. شبیه‌سازی مدل در ماه‌های ژانویه, فوریه. مارس و آوریل از تطابق خوبی برخوردار 
بوده و جریان شبیه‌سازی شده به‌خوبی جریان مشاهده شده را دنبال می‌کند. مقدار حجم رواناب محاسبه شده 


۵ میلیون مترمکعب و متوسط روانات محاسبه شده ۰۸:۹ مترمکعب در ثانبه می‌باشد. مقادیر دبی و حجم 


طمتاحصتصصجهاع ۶ اصمتم ]0۵1 16 1 
6 ماه ۷ 2 
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۱ جغرافیا و مخاطرات محیطی شمارة اول 


رواناب مشاهده شده و محاسبه شده بسیار به یکدیگر نزدیک هستند. پس از آن اجرای مدل با فاکتور درجه- 
روز کلاسیک و باز هم استفاده از دو ضریب فروکش فصلی در شبیه‌سازی موفق‌تر بوده و دارای ضریب تبیین 
۹:9 برابر با ۰/0۲ و درصد اختلاف حجمی (,) برابر با ۱7/64 بود. این روش دارای کم‌برآورد در مقدار 
جریان شبیه‌سازی شده در ماه‌های نوامب دسامبر و فوریه و مارس بوده است. مدل ۹1۷ ضعیف‌ترین 
عملکرد را در اجرای مدل با استفاده از فاکتور درجه-روز کلاسیک و یک ضریب فروکش سللانه برای کل 
دوره شبیه‌سازی را از خود نشان داد در این اجرا ضریب تبیین (*18) و درصد احتلاف حجمی (10) به‌ترتیب 
پرابر با ۰/6۳ و ۱۹/۲۲ بودند. در این روش مدل در ابتدای ماه نوامبر و ماه‌های پایانی دوره شبیه‌سازی دارای 


کم‌برآوردی می‌باشد و استفاده از این روش در شبیه‌سازی رواناب ذوب برف توصیه نمی‌شود. 
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3 1.5 
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رواناب محاسبه شده با فاکتور »جدید و دو ضریب فروکش فصلی نت 
رواناب محاسبه شده با فاکتور کلاسیک و دو ضریب فروکش فصلی مس . . سب 
رواناب محاسبه شده با فاکتور جدید و یک ضریب فروکش سالانه - - - - 
رواناب محاسبه شده با فاکتور »کلاسیک و یک ضریب فروکش سالانه مموو وود ممممووه 


شکل ۶- رواناب اندازه‌گیری و شبیه‌سازی شده توسط مدل 518۷1 برای سال ۲۰۱۱۲۰۱۲ در حوضه کارده 


جدول *- مقادیر معیارهای ارزیابی دقت مدل 1811 در اجراهای مختلف در سال ۲۰۱۱۲۰۱۲ 


درصد اختلاف 
حجمی 1 ضریب تبیین "1۴ روش‌های اجرای مدل 
1۷ ۰/۷۲ رواناب محاسبه شده با فاکتور 0 جدید و دو ضریب فروکش فصلی 
۱/۹۹ ۰/۳ رواناب محاسبه شده با فاکتور 0 کلاسیک و دو ضریب فروکش فصلی 
۷/۷ ۰/۱ وتاب سخاسبه نشده با فاکتور 09 ید .ی یک ضیریب فز وکین لاه 
۱۹/۳۲ ۳ دواتات مطانسه شاه با قاکتوی و6 کازسک ور یک ربب قر رگن مبالانه 


درستی سنجی مدل: شکل ۵ هیدر وگراف‌های روانابت اندازه گیری شده و شبیه‌سازی شده توسط مدل 
در حوضه آبریز کارده در سال ۲۰۱۳ ۲۰۱۲ که جهت صحت‌سنجی مدل اجرا شده را نشان می‌دهد. 


همان‌طور که در (شکل ۲ نیا مشاهده می‌شود هر این سال شاهد کاهش محسوس سطح پوشش برف و در 


سال سیزدهم مدل‌سازی رواناب ذوب برف توسط مدل 71 و برآورد پارامتر فاکتور درجه-روز با استفاده از ۱۵ 


نتیجه کاهش حجم آورد و نوسانات دبی مشاهده شده نسبت به دوره واسنجی هستیم که خود باعث کاهش 
طقت تسا یرتاب کوب بر فرط ال کرو اباب امین سارهاش ارزیانی خر هر سل 
۴ مورا 4۴ریت یت 0و درضت اکباای ی( سا اسای بو با برش 
محاسبه رواناب با فاکتور درجه-روز جدید و استفاده از یک ضریب فروکش سالانه به ترتیب برابر ۰/۵۱ و 
۸ بوده و دارای بهترین عملکرد در شبیه‌سازی رواناب توسط مدل 91۷ از میان روش‌های موردبررسی 
در این تحقیق است و نتایج حاصل از دوره واسنجی را تأیید می‌کند و تطابق خوبی را برای دبی‌های اندازه- 
گیری شده و شبیه‌سازی شده در حوضه کارده نشان می‌دهد. مدل در ماه‌های ابتدایی شبیه‌سازی, نوامبر تا 
نیمه دسامبر و ماه‌های پایانی شبیه‌سازی مارس و آوریل دارای کم‌برآورد بوده است و مقدار حجم رواناب 
همقل 6۳۱ میرن یکب و سط رالات اقلاوه کر ی ۱/۳6 ری کف بر کته ید 
باشد و مقادیر حجم رواناب محاسبه شده و متوسط رواناب محاسبه شده به‌ترتیب ۵/14۷ میلیون مترمکعب و 
۲۷۱ مترمکعب بر ثانیه است. پس از آن نیز اجرای مدل با فاکتور درجه-روز جدید و دو ضریب فروکش 
فصلی در شبیه‌سازی موفق‌تر بوده است. این اجرا دارای بیش‌برآورد در ماه‌های ژانویه. فوریه و آوریل می- 
تاکز ریب یی ارو فرصت ات یی 0ب کب را 9 ۳۲ انس سرا 
مدل با استفاده از روش محاسبه فاکتور درجه-روز کلاسیک و به‌کارگیری دو ضریب فروکش فصلی در طول 
روند شبیه‌سازی ضعیف‌ترین عملکرد مدل را از میان روش‌های به‌کار برده شده در این تحقیق به نمایش می- 
گذارد. در این اجرای مدل ضریب تبیین (87) برابر با ۰/۳۳ و درصد اختلاف حجمی (,2) برابر ۱۵/۳۲ می- 
باشد. در این شبیه‌سازی مدل دارای کم‌برآورد در ماه‌های نوامبر و دسامبر می‌باشد و نتایج حاصل از دوره 


زمان (روز) 0 
۰/۳3 ۰ 03/01/2013 02/01/2013 ۰ 01/01/2013 اناب12/451:/2001 شد2 11/04/2۵ 
رواناب محاسبه شده با فاکتور »جدید و دو ضریب فروکش فصلی دجوجد2د 
رواناب محاسبه شده با فاکتور ,کلاسیک و دو ضریب فروکش فصلی س. . . سب 
رواناب محاسبه شده با فاکتور جدید و یک ضریب فروکش سالانه - - - - 


شکل ۵- رواناب اندازه‌گیری و شبیه‌سازی شده توسط مدل 518۷1 برای سال ۲۰۱۱۲۰۱۲ در حوضه کارده 


۱۹ جغرافیا و مخاطرات محیطی شماره اول 


جدول ۷- مقادیر معیارهای ارزیابی دقت مدل 81811 در اجراهای مختلف در سال ۲۰۱۳ ۲۰۱۲ 


درصد اختلاف حجمی ضریب تبیین 


:10 9 روش‌های اجرای مدل 

۱/۳۷ ۷ رواناب محاسبه شده با فاکتور 0 جدید و دو ضریب فروکش فصلی 
۱۹/۳۳ ۳۳ رواناب محاسبه شده با فاکتور 0 کلاسیک و دو ضریب فروکش فصلی 
1۳۸ ۰/۵۱ رواناب محاسبه شده با فاکتور 0 جدید و یک ضریب فروکش سالانه 
۱/۲ ۳۷ رواناب محاسبه شده با فاکتور 0 کلاسیک و یک ضریب فروکش سالانه 


کشور ایران در منطقه گرم و خشک قرار دارد و میزان حجم توده‌های برفی در حوضه‌های آبریز اين گونه 
مناطق مانند حوضه آبریز کارده محدود می‌باشد و اين یکی از دلایل کاهش دقت شبیه‌سازی توسط مدل 
است. میزان متوسط بارش برف در حوضه آبریز کارده در سال ۲۰۱۲۲۰۱۳ برابر با ۳۰/۰۵ سانتی‌متر بوده و 
از دلایل کاهش دقت مدل در دوره درستی‌سنجی نسبت به دوره واسنجی می‌توان به کاهش ۲۷/۰۲ درصدی 
بارش برف و کاهش آورد سالانه اشاره کرد. همچنین خطای انسانی در اندازه‌گیری رواناب را از دلایل عمده 
آن دانست. چنانکه قربانی‌زاده و همکاران (۱۳۸۹) نیز در مطالعه خود که بر روی رواناب حاصل از ذوب 
برف با کمک مدل 518 در حوضه آبریز کارون انجام دادند و دی و همکاران" (۱۹۸۹) در شبیه‌سازی ذوب 
برف در پاکستان به عدم تطابق دبی‌های حداکثر مشاهده شده و شبیه‌سازی شده اشاره کردند و بیان کردند که 
این ضعف ناشی از اشتبامات در داده‌های مشاهداتی و بی‌کیفیتی آن‌ها می‌باشد. نتایج به‌دست آمده از شبیه- 
سازی رواناب حاصل از ذوب برف در این مطالعه با نتایج حاصل از تحقیقات میریعقوب‌زاده و همکاران 
(۰۱۳۹۰ اربابی سبزواری و نجاری (۱۳۹7) و آرتیمانی و همکاران (۱۳۹۸) مطابقت دارد. همچنین مطالعه 
قاسمی و همکاران (۱۳۹۲) که نوسانات رواناب حاصل از ذوب برف تحت تأثیر پدیده تغییر اقلیم در دهه- 
های آینده را موردبررسی قرار دادند و در واسنجی مدل ضریب تبیین و تفاضل حجمی را به‌ترتیب ۰/۷۱ و 
۵ درصد به‌دست آوردند نتایج به‌دست آمده را تأیید می‌کند. آن‌ها نتایج حاصل از شبیه‌سازی واسنجی 
مدل را قابل قبول ارزیابی کردند. نتایج مطالعه لی و ویلیامز" (۲۰۰۸) بر روی رواناب حاصل از ذوب در 
حوضه بدون آمار تاریم در چین که مقدار ضریب تبیین را در دوره واسنجی ۰/۶ و برای دوره درستی‌سنجی 
۸ و ۰/۵۱ به‌دست آورده‌اند نیز با نتایج تحقیق حاضر هم‌خوانی دارد. نتایم حاصل از محاسبه فاکتور 
درجه-روز با استفاده از روش جدید (معادله شماره ۶) و روش محاسبه کلاسیک (معادله شماره ۱) و مقایسه 
نتایج آنها نیز با مطالعات مارتینک (۱۹۸۹)» کوستاس و رانگو (۱۹۹۶) و بروبکر و همکاران (۱۹۹7) وفاخواه 
و همکاران (۲۰۱۶). زیمنس و همکاران (۲۰۲۱). دلاور و همکاران (۱۳۹۰). ابراهیمی و همکاران (۱۳۹۵). 


۰ 16۷ 1 
محصهتا از ۷۷ 0صح 7[ 2 
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رشیدی و همکاران (۱۳۹7). مطابقت دارد. نتایجم حاصل از اجرای مدل با روش‌های مختلف برآورد فاکتور 
درجه-روز در دوره‌های واسنجی و درستی‌سنجی حاکی از برتری روش جدید برآورد فاکتور درجه-روز 
(معادله شماره ۶) در شبیه‌سازی رواناب ناشی از ذوب برف نسبت به روش کلاسیک محاسبه آن (معادله 
شماره ۱) است. نتایج به‌دست آمده نشان می‌دهد با استفاده از روش جدید محاسبه فاکتور درجه-روز و 
محصول تابش خالص هشت روزه (۴۲۲۳۱/۲-۴(-5۳5)) ۱520 می‌توان با روشی ساده‌تر و کم‌هزینه به 
مقادیر قابل قبولی از فاکتور درجه‌-روز (0) دست پیدا کرد. همچنین نتایج تحقیق حاضر نشان می‌دهد 
ضریب فروکش جریان (1) پارامتری حساس در شبیه‌سازی رواناب ذوب برف توسط مدل ٩1۷۲‏ است و 
تأثیر به‌سزایی در تطابق جریان اندازه‌گیری شده و شبیه‌سازی شده دارد و در نتیجه تعیین دقیق آن می‌تواند 
نتایج بهتری با توجه به سنجه‌های کارایی مدل به‌دست آورد. نتایم حاصل از تحقیقات وفاخواه و همکاران 
(۲۰۱۶) نیز اين مطلب را تأیید می‌کند. نتایجم حاصل نشان می‌دهد با توجه به اينکه مقدار ضریب فروکش 
جریان متغیر بوده. در دوره‌هایی که نوسانات جریان و آورد رودخانه زیاد است مدل 830۷1 می‌تواند با استفاده 
از دو مقدار ضریب فروکش فصلی در شبیه‌سازی نقاط اوج موفق‌تر عمل کند و در دوره‌های که جریان 
یکنواخت بوده و نوسانات آن کم است. استفاده از یک ضریب فروکش سالانه عملکرد بهتری دارد. 


۵- نتیجه گیری 


با توجه به نتایج حاصل از اين مطالعه و اجرای مدل با روش‌های مختلف در دوره‌های واسنجی و 
درستی‌سنجی. بهترین روش مدل‌سازی جریان در دوره واسنجی» روش محاسبه فاکتور درجه-روز جدید 
(معادله شماره ۶) و کاربرد دو مقدار ضریب فروکش جریان فصلی با ضریب تبیین ۰/۷۲ و درصد اختلاف 
حجمی ۶/۱۷ بود. در دوره درستی‌سنجی نیز روش شبیه‌سازی رواناب با فاکتور درجه-روز جدید (معادله 
شماره ۶) و یک ضریب فروکش سالانه بهترین نتایج را از نظر معیارهای ارزیابی دقت مدل ارائه داد. ضریب 
تبیین در آن ۰/۵۱ و درصد اختلاف حجمی 1/۳۸ بود. ازاین‌رو می‌توان نتیجه گرفت محاسبه مقادیر فاکتور 
درجه-روز (0) با استفاده از معادله جدید ابداع شده در اين مقاله (معادله شماره ۶) و تعریف مقادیر آن برای 
مدل 581 به مدل کمک می‌کند تا عملکرد بهتری نسبت به شیوه کلاسیک محاسبه فاکتور (0) با استفاده از 
داده‌های چگالی برف اندازه‌گیری شده در حوضه آبریز (معادله شماره ۱) ارائه دهد. ترکیب مولفه تابش در 
معادله محاسبه فاکتور (0) نیز نتایج مطالعات گذشته مانند کوستاس و رانگو (۱۹۹۶) و بروبکر و همکاران 
(۱۹۹7) در ترکیب مولفه تابش در محاسبه عمق ذوب برف در مدل 580۷ را تأیید می‌کند. محصول تابش 
خالص استفاده شده در این مطالعه برای محاسبه فاکتور (0) به روش جدید -2۲۳1/۶(-05۳85۲5) 80 


(ظ نیز از دقت بسیار خوبی برخوردار بوده است و می‌تواند به نحو قابل قبولی برای شیبه‌سازی رواناب 


۱/۸ جغرافیا و مخاطرات محیطی شمارة اول 


حاصل از ذوب برف در حوضه‌هایی که با کمبود داده مواجه هستند مورد استفاده قرار گیرد. همچنین نتایج 
معیارهای ارزیابی دقت ارائه شده نشان می‌دهد ضریب فروکش جریان نقشی اساسی در تطابق جریان‌های 
اندازه‌گیری شده و شبیه‌سازی شده ایفا می‌کند و همانطور که در مطالعات گذشته مانند وفاخواه و همکاران 
(۲۰۱۶) آمده پارامتری حساس و موثر در شبیه‌سازی رواناب ذوب برف توسط مدل 518۷1 است. نتایج 
تعاماز از آنخ قوعش شان داه ایج شرب مر اند صورت الا با فان راخ ان 9۳01 مرت گزود 
و تعیین دقیق آن در بهبود عملکرد مدل کمک شایانی می‌کند. ضریب فروکش جریان مقداری ثابت نیست و 
با کاهش دبی مقدار آن افزایش می‌یابد و بالعکس. شبیه‌سازی رواناب ناشی از ذوب برف در حوضه آبریز 
کارده نشان می‌دهد برف قسمت قابل توجهی از آورد رودخانه در فصل زمستان و بهار را تشکیل می‌دهد. 
ذوب برف از اواخر ماه فوریه در حوضه افزایش می‌یابد و می‌تواند در صورت همراه شدن با رواناب حاصل 
از بارش باران منجر به وقوع سیل به حصوص در اواخر ماه مارس و ماه آوریل گردد که نیازمند برنامه‌ریزی و 
مدیریت کارآمد در جلوگیری از خسارات ناشی از سیلاب‌ها و استفاده بهینه از اين منابع آبی می‌باشد. با 
توجه به نتایج اين پژوهش مطالعه بر روی رواناب ناشی از ذوب برف در حوضه‌های آبریز امری ضروری 
است و می‌تواند در مدیریت بحران و برنامه‌ریزی برای مقابله با مخاطرات و خسارت‌های حاصل از سیلاب- 
ها اف راز دوب هدور ازا ضا رفستاق ز ابایل قصا هار نک تایه ای وی بایان 
کرده و از بروز خسارات جبران‌ناپذیر جلوگیری به عمل آورد. اين‌گونه مطالعات بر روی زمان‌بندی و میزان 
رواناب حاصل از ذوب برف به مدیران در برنامه‌ریزی و مدیریت کارآمد برای تأمین و توزیع آب مورد نیاز 
بخش کشاورزی در فصل‌های گرم سال. جهت جلوگیری از مخاطرات ناشی از کمبود آب به‌خصوص در 
سال‌هایی که کم‌بارش بوده پا به علت گرمای زودهنگام هوا ذوب زودهنگام توده‌های برف اتفاق می‌افتد. 
کمک شایانی می‌نماید. با توجه به اينکه کشور ایران در منطقه خشک قرار گرفته مطالعات مستمر این‌گونه 


منابع برای مدیریت موّثر و کارآمد و جلوگیری از هدررفت آب ضروری می‌باشد. 
این مقاله از رساله دکتری با عنوان ارزیابی اثرات تغییر اقلیم بر توده‌های برف و رواناب سطحی (مطالعه 
موردی: حوضه آبریز کارده) در دانشگاه فردوسی مشهد استخراج شده است. 
کتاینامه 
ب ارتیمانی» محمد مهد ی؛؟ زینی‌وند. حسین؛ طهماسبی‌پور: ناصر؛ ۳۹/۸ ارزیابی کارآیی مدل ۷ و 11۳9۷ در شبه - 


سازی رواناب ناشی از ذوب برف در حوزه آبخیز بوجین. تحفیقات مناب عآب ایران» ۵( ۲۶۱-۲۲۸. 
2 ۵ ۱)0://۷۷۷۷۷۷ 
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- ابراهیمی» راضیه؛ حمزه. سعید؛ معروفی. صفر؛ ۱۳۹۵. مدل‌سازی سطح پوشش و رواناب ناشی از ذوب برف با 
استفاده تلفیقی از مدل هیدرولوژیکی 5101 و تصاویر ماهواره‌ای. مهندس یآبیاری ‏ وآب ایران ۱ (۲۳ ۷۷-۷۲ 
0 و( 
- اربابی سبزواری, ازاده؛ نجاری مهربان. فرنیاز؛ ۱۳۹۹. نقش مدیریت منابع اب در توسعه کالبدی با استفاده ازمدل 
مطالعه موردی حوضه آبخیز طالقان. برنامه‌ریزی و توسعه کالبدی, 1(۲ (پیاپی 6۸ ۸۵-۷۸ 
یه ۰۱۰/۸ ۱۱۱۹۱ 
- جهانبخش اصل. سعید؛ دین‌پژوه» یعقوب؛ عالی‌نژاه محمدحسین؛ ولی‌زاده کامران. خلیل؛ پرهیزکار مرتضی؛ ۱۳۹۵. 
شیبه‌سازی رواناب ذوب برف در حوضه شهرچای با استفاده از مدل 5181 جغرافیا و برنامه ریزی محیطی (مجله 
پژوهشی علوم انسانی دانشگاه اصفهان» ۲(۲۷ (پیاپی 1۳)): ۱۶-۱ 
2 09/۷۷ 
- دلاون مجید؛ مرید. سعید؛ نیکبخت. نجمه؛ ۱۳۹۰. شبیه‌سازی توزیعی ذوب برف در حوضه‌های کوهستانی فاقد داده 
(مطالعه موردی: حوضه امامزده داوود). تحمیقات ماب عمآب ایران ۷( (مسلسل ۲۲)): ۵۰-۶۱. 
:2 و 009/۷۷۷ 
- رشیدی. محمد؛ حاجی بیگلو» محبوبه؛ سربازی» محبوبه؛ قادری» محمد؛ ۱۳۹۲. برآورد رواناب حاصل از ذوب برف 
در حوضه‌های آبریز استان خراسان شمالی با استفاده از مدل ۷۷:0۹[2[۷ (مطالعه موردی حوضه آبریز دربند 
سملقان). علوم مهندسی ‏ وآبیاری (مجله علمی کشاورزی) ۲(۶۰). ۱۷۱-۱۵۹ 
2 3 ۷۷۷۷۷۷//: ۱0 
- طالبی اسفندرانی» سمیه؛ علوی‌پناه. کاظم؛ علی‌محمدی سراب. عباس؛ روستاء حسن؛ ۱۳۹۰. جداسازی ابر از برف 
در تصاویر ۷10115 با استفاده از الگوریتم مق ٩20‏ و الگوریتم ماسک ابر. مجله سنجث رازدور و 5آت) ایران. 
۳ ۷۱-۹۰ 
8 ۷۷۷۷//: 05 
- عالی‌نژاد. محمدحسین؛ جهانبخش. سعید؛ خورشیددوست. علی‌محمد؛ ۱۶۰۰. شبیه‌سازی رواناب ناشی از ذوب 
برف حوضه قره‌سو با مدل ]95[۷. جغرافیا و برنامه‌ریزی. ۰۷۱(۲۵ ۷-۸۷ 
58 ۵۱۱۱۵۵۵0 ۱605۰/۱۷۷ 
- قربانی‌زاده خرازی. حسین؛ صدقی. حسین؛ ثقفیان؛ بهرام؛ پرهمت. جهانگیر؛ ۱۳۸۹. پیش‌بینی توزیع زمانی جریان 
رواناب ناشی از ذوب برف در نیم قرن آینده تحت شرایط تغییر اقلیم. مهندس یآب, ۱(۱ ۵۲-4۳. 
2 1( 
- قاسمی. الهه؛ فتاحیء ابراهیم؛ بابائی ام‌السلمه؛ ۲ بررسی نوسانات رواناب حاصل از ذوب برف تحت تأثیر 
پدیده تغییر اقلیم در دهه‌های آینده. نشریه پژوهش‌های اقلیم شناسی, ۱۳(۶ ۰۱۱۲-۱۲۲ 
لصط. 14153 _مامنا1۳/2۲. ما1 ما ۱۷۱۱۷۱۸//:ح۸ظ 
- قنبرپون محمدرضا؛ محسنی ساروی» محسن؛ ثقفیان. بهرام؛ احمدی. حسن؛ عباسپون کریم؛ ۱۳۸۶. تعیین مناطق 
موثر در انباشت و ماندگاری سطح پوشش برف و سهم ذوب برف در رواناب. مجله منابع طبیعی ایران. ۳(۵۸) 
۵۱۲-۵. 6 ۰۵22 1005۰/۷۷۷( 


۳۰ جغرافیا و مخاطرات محیطی شمارة اول 


- قنبرپور میرحسن؛ قنبرپور» محمدرضا؛ حبیب‌نژاد روشن. محمود؛ ۱۳۹۰. مدل‌سازی جریان ناشی ذوب برف با 

استفاده از مدل هیدرولوژیکی رواناب حاصل از ذوب برف (مطالعه موردی: حوضه آبخیز سد کرج). تحیقات 
مناب مآب ایران, ۳(۷ ۵۲-۶۰ 

۵۵5۰/۱۷۱۷ ۰۱۳۸۱۲۲۰۱۲/۵۲)۵1۵_ 16219. 

- میریعقوب‌زاده» میرحسین؛ قنبرپور محمدرضا؛ ۱۳۸۹. بررسی کاربرد نقشه‌های پوشش برفی حاصل از تصاویر 

ماهواره‌ای ۷101915 در مدل‌سازی رواناب ذوب برف (مطالعه موردی حوضه آبریز سد کرج). مجله علوم زمین, 

29 ://0(.۵/10.22071/ 2 2 .۱۶ ۱-۱۶۸ ۹ 


- نجف‌زاده. رضا؛ ابریشم‌چی. احمد؛ تجریشی. مسعود؛ طاهری وه حمید؛ ۰.۱۳۸۳ شبیه‌سازی جریان رودخانه 
با مدل ذوب برف (638[۷). مجله آب و فاضلاب ۱۵( ۱۱-۲. 


لصاح 2487 مامتاه/1 ۰۸۸0۵۱۵ ۸05۰//۱۷۱۷۱۱۸ظ 


ج معصقطه ماقصتاه ۵۶ ]قمحا ,92016 رفظ مفطوطصمطهل ی ولا مطمطفهم‌مانا رب ۷۰ بله[م۸118 - 

۰ 6288) ۵0615 60۷ ۵۶ باصتهاهعصن معط عصنامممع مفصا عصت‌له) وا امه صم ۳۵0۴ 

,(5)1 ,۹0100665 هم هه اعدا ۵۶ افصتول مصتلمون) مهممممب ,رفتصصرتا مد صتفقظ رفطم‌طهط 
۰ .0161206 60۶۲00۵6210-5 //:۸۵۹ظ 200-211۰ 

۲ مصو فطع مفصا قتامرصا مهتم ممتاه10601001 .۷۷۲۰,1996 ,قه)وناکا عک رنه ,عطق ویک رتعلهانارظ - 
 60۵6://001.0۲8/10.1002/)5101(1099-‏ .1329-1343 ب,(10010 ,۳۲۲۵68685 ۳۱۷۵۲۵1 ,061مظ. ما۲ 
۳464<3.0.00:2-۷۰ ۷ ]199610(10:10>1329::۸]1-۳) 1085 

۶۲ 210۳ 2 ۲0۲ [00م گمصیاتآممه مد 0 ادع) ۸ .1989,.ظ ,معصقک ‏ ریک ۷۰ مقحطتقطاه ر.ظ رلامر[ - 
.167-۰ «(2063 ,16562100 زهروازی ۱۶۱۵۱ ۰ و ۱ 11 929 
۰( 00( 

ملفصتاه عصاووعوعوظ ,۱2019 رطاتصه عک بیط رل مقلل۱۷0 بب بن) بعاععاو بیظ ,21020 بلط بط ,و۲112 - 
ملک عمو۱۲] عطا ۵۶ مصتعهها. همهم امصومصی که دنه ماهته وه مامهم‌طظا ععصقطه 
۰ 16/[.0[10.20153 10.100 /05://001.0۲82 .325-546 ر3 ,۲۱۷۲0۱02۷ ۵۶ امدرتان0 صنفجها عصهری) 

2 1 امهمصص ۶مصیت تافص بس«مصو ۵ ممتاجعتانموظ .۷۵2016 مصفلطعم0هه عک بیط بل2تا۳1۳0 - 
74-1۰ ,(4)2 ,۲۱۵۲۵۵0108 )م۱1۲0 0ج ممممزمو06۵) ۵۲ اقصتامل مصق! ۵ صتعجها تامصتع)صتا۵ظ 
۰ 1/0( 

2001162010 هه اممصصع۷ب۲0مرح .2019وناط بلنا ورگ ر9018 ویک معطهلا ویک بان ول بل بت) ,۳120 - 
36798 ,26 ,۴۵6۵۲0 ممتناا۳۲0 20 عمممهی لم)ممصصمرم1صظ .عممنععع مصتمرله ما 9 عط ۵۶ 
.6۵5://001:0182/10.10017/811356-019-06814-۰ :36811 

1۰ مصن )فصو مط ز10 اعع۵0 عم عامرصصلو ۸ 1994 ,0عصق ‏ ر.ظ ۷۷۰ رفهادتاک - 
۷۰ ,1515-15327 ,(30)5 ,5۵56276 وعم‌عنا0ومع تمه ۱۷۷ 

,۷۵۲۵۲۹0۵0 صتمتصنامصه نع صرح ما عصتالم0مصط مصند اممهبمصرو ۷۷۰2008 ۱۷۲۰ رفصصهتا۷۷1 ی ربعر بان - 
3931-0۰ ,(22)19 ,0۲06۵6۵8565 ۳۱۷0۲01021021 .قصنط رصتفحظ صصتته 1 

0۹://01.012/10 1 (۰ 

بصع و موه (8۳۴) اع0محط من امجصچمصه ,2008 رعاته‌حام؟ ک ریق ر0عصق؟ رل ر6ع۷]2۲110] - 
,6003 ۱۷ ,۱0665 وق راصنا متهاگ م۵ م۱ رمعصفگ آمتممصتهم۲ ۵۲8۵202[ باون 
060۴۵0100.۰ ۳/۷۲5 20/9 ۷ 00ع. تافص وتا //:05ظ :وتا 

اجه مه عصتنل مان عامصصتوت۱ مضه ماهر ]امه مصی ‏ تنامط‌م)ا ۲]0‏ .1.1989 رمعمتایم]۷ .- 
عمصتلقظ عمط ۵۶ مهن رقصمتامانه ۷‏ علمهاع .مه ومع . ۹۵ مه طمتاقآهاه 
,19-8۰ ,183 اطن۲ فها ,معصقط ‏ به. بها . 0هال0ع . ,1989 ۷۲2۲ . رلصها ها . رصتاتدمم‌ص 5 
06016۰ 2۲6۲/۲۷/۲۹ ۷ /۵00: افرص ومامام //۵9۰ 

۷۲ص ما الازصاهته۷ لهز۵محصمامل5۵2 .2020بخ روهظ ی ره ره 070۵20صقظ .۷ متصهططمصرن لا - 
1۳1(۷(۰) ۱۵061 من )امه مصی 0ماماتتولل ت«الفتاهمو موه ماد و0۲ صم وتعاعصهتهم 


سال سیزدهم مدل‌سازی رواناب ذوب برف توسط مدل 71 و برآورد پارامتر فاکتور درجه -روز با استفاده ازع ۳۱ 


عامصصع ۵۶ جامزهمو صقل4ظ1 مطا ۵ لمصنمل ,فرظ لمطمهصننم رمتفهط عصوه‌طارنا ما فتاه عقوم 2 
.10007/512524-020-01215-۰ 8/10۰ ۳۸۵://:۰001۰01 :325-340 49 رعطلفممو 

۶ موجه صح مملاهعتآص۸ .۲.,2014 ,۱۷۲ بظ10۷7ظ وی یک ,۳۲۵ وه ۷۷1112108 ویک .ظ ,۳2002۷ - 
۷ مطلودا ۳۱۱12۵ صعاعجء معط ما متفجها بعن عمصصح 1 فط مر اممصهه مصتار )م90 ۵ 
«(28)21 ,۶۵6۵85 ۲۱۷۵۲0۱۵۵162 .موجه مم)هاتصتفوه فاهل ۷0۵ ماتقن ۷۲۵۵۸8 صتقطل) 
۰ 10۰1002/۳7 /۳۸/05://001.0۲2 5337-5353۰ 

مفصتوله جح مظ ممتامع[0۳0 ۲مصنن؟ :1.2021 ۲۲۵۱۷۹۵ عک ریک بک مقطاوه‌تط رل مقعلطاانا ریک رفصممر - 
۰ ,(13)9 ,۷۷۵6۲ .۵061 مصن )ام مرو ه ۵۶ مماه)ناممه تعلهممت تفع ما 60وم2 ۷۱ 
۰ و( 

2 1 وصتاام0مصظ مصن اممووجمصو. .2015ریه ره وطل070۵۵ص2ظ ی ره ب2072ظ و.؟ به0ام8602 - 
۰ 629122 عامصصمک مممتععر مه رملمصنا صمافهه ۵۶ پمعصامتهه عصمممت ظ 902168 
3۰ 00۵5://001.0۲8/10.1016/1.15286.201 20-35۰ و1 ,اصعحصمما حظ مه 50060 

۰ ,۲161۳۷0199 ی بظ ,0عصقط ویظ مفقلاط بلا مففصطهل رل ,112168240 بب ) ,906010 - 
ملک معط ۵۶ ماو موه مص اامصومصی ممتاام0مصه رم فامالمن مرو فانانم عمتاهتاه۷ظ 
,0۲۱۵08ص01عع مه ممتاهموهان وه نامه گم لحصنمر امممتامممام .عماه 620 مهن 
۰ .160/1 00://0:001.012/10.10 .234-243 ,63 

ملجع۵0 ۵ صمتاهعلانظ 2014بیظ .1 بهو۷]0 تک ب.ظ .ظ بلفحططظ ویک م2021 وب متطفطفطه - 
۶مصن امصعمصو . مصتاملمصهر. صا دوصقم مهن مگ مامصمه۵1ا۵۵0د . عصلعقصظا ممتاتاا0وع۳ 
1-۰ 8 ,عصلقصع0 عامصصفک معناممخ ۵۶ حدم مصقا رصتفدها._ هم تقتاصوع فط صا ووععمم 
۰ ».7/1 ۷۷://0::01.0182/10.111 

عصلما ممتامصتای کمصیاتتامموهمصو ,2014ییک ره رطقمهمد ها عک بیظ ماما ۷ مطفقط2121 ۷ - 
993-۰ ,73 مومعمصماهی طاتقظ مماص بمطفماهس مقطومله 1 ص اممصه 9۳۴ صمتاهل1۵0 
-۷۸0://0::01۰0۲8/10.10017/512665-014-3449-5 

6 ۸۱۱۵۲ عصلودا ممتلح‌ممهع گممنار مه مم۵ 9۵۵0 1 ۲.۷],,201 رد1260 عک بت ب۱0۷0۷10 ۱۷ - 
۰ )1۵۱۵0 .صه)فتممطعگه رلعطوع2 ۱۷۷ 0محصاقظ بممهن۲ا فط ما فممتاه و۵۵ ۸1۵۲0۵۷۵۷۵ 
۰ 1 ۳05://::001.018/10.1016/1 .3313-3321۰ ,115 ,۷1۲08 

۵ .۳۲۱.,,2018) بلار ‏ ما2 بلنا .۷۷ ,206082 کر بعظ۲۵ میک م202108 وک بتا200 رل بتالل ملگ ر2410 - 
عطا ۶ تامهم اصمله۷ ما عصتاهطاللمه بم صط)تمعله مامتان لماوع 0۳۵2۲6۵891۷۵ 
,470-۰ . ,566 ۳۱۷۲0۱02 ۵ .. لقصتامژ. .قمع مت امصظجمصو 
۰ .10۰1016 /۱)0۹۰//001.0۲2 


